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XI. On certain Definite Integrals occurring in Spherical Harmonic Analysis and on 
the Expansion^ in Series^ of the Potentials of the Ellipsoid and the Ellipse, 

By W. D. NiVEN, M.A., Felloiv of Trinity College, Camhridge.. 

Coramunicated hy J. W. L. Glaisher,, M.A., FM,S., 



Beceived April 3j—E-ead May 1, 1879. 

§ 1. The object of tHs paper is to explain a method of dealing with a class of results, 
some of which are of frequent occurrence and some of considerable importance, in the 
solution of Laplace's equation in series. 

The fundamental theorem on which the method depends is expressed by the inte- 
gration of 

over a sphere whose centre is the origin of coordinates and whose radius is R. This 
sphere we shall call the sphere of reference. 

By change of axes the aboye integral takes the form 



27rR e'^^^'-'^'-'y'^dx- 



and its value is 



27rR 



gBV(a^+/3=^+y=^)__g-EV(a=^+^2^.y2). 



9 • • » • •■ » \ -^ / 



or, in series, 

l+|V+^'+/)+ . . . +2^^ . . (2) 



47rll^^ 



It follows that if V be any function whose value may be expressed for all points 
within the sphere by a convergent series, according to Tayloh's theorem, or by the 
symbolical form 



^ dx dy dz ^ 



then the integral ||V(iS taken over the surface of the sphere is 



} 



^^^12T;TiU+c^,+ jJ^o (3) 



3 c 2 
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The corresponding theorem in two dimensions, which we will also employ, is this :— 
The integral 

taken round the periphery of a circle of reference, is 

CO J^3I 

27rRS;T^r77T,(aH/3^)' .....«*.. ^ (4) 

Hence it follows that 

taken round the same periphery, is 

The theorems expressed by (3) and (5) are employed below in a variety of cases con- 
nected with spherical harmonics. It is therefore necessary to specify at the outset 
the notation to be used. 

We denote, according to custom^ the colatitude and longitude of a point referred to 
the axes of reference by and <^ : cos ^ and sin Ohj \l and v respectively. We denote 
the Legendre's coefficient or zonal harmonic of the i^^' degree by the symbol P,- so 
that P^ is defined by the equation 

6} 1 

V I XT ^' " ' I """""' i- / I -J » • * «■ « " » » « » \/ 

According to Clebk Maxwell's theory of poles the general expression for the har- 
monic of the i^^' degree is given by the equation 

.'! y.^/^nv+i— ^^ (7) 

where — means differentiation with regard to an axis L 

dh 

In the harmonic (7) the i poles are the points where the i axes of differentiation cut 
the sphere of reference. In the harmonic (6) the f poles all coincide with the point 
where the axis of z cuts the sphere. If the i poles coincide at any other point the 
harmonic will be denoted by Q^-, 

In the tesseral and sectorial system we have ventured to depart from the usual 
notations by denoting them as follows :— 

ii (,»=(-i)v«(|+|)i;'^ ... ... (8) 



n 



i, a-j^^i-lp" 3\J-^~^j^ ; ....... (9) 
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where y= a/ — 1? ^=^+^2/? y)=^—'jy* It will be observed that (8) is not in accordance 
with the general definition (7), for if we put <t=0 we get (^, o)=:2P^-. 

As regards the name of Spherical Harmonic introduced by Thomson and Tait, it 
seems desirable that the names of Legendre and Laplace should be in some way 
retained in connexion with the different kinds of harmonics with which they are asso- 
ciated. From the point of view of the former writers there is perhaps some difiiculty 
in doing this. We shall however continue to give the name of Laplace's coefficient 
to what Thomson and Tait call the Biaxal Harmonic or the formula which expresses Q^ 
in terms of the harmonics referred to fixed axes of reference. With the definitions 
given above, this is 

Qi=V^/+, &c., + .^^fjj (i, a){i a)'+\i, cr]\i, crj}+ > • (10) 



in which any accented function is derived from the corresponding unaccented function 
by the substitution of the coordinates of the pole of Q^ in place of the running co- 
ordinates. 

The following are, in the first instance, the integrals to which the method of this 
paper will be applied : — 

(1) f^P.c^S 



(2) j'f(l-^^)-p,,c?S 

(3) |f(l-/x^)-^P,,c^S 



(4) 



X. /nX. qJL f\A-Kj 



(5) {[(p, a)(^, P){r, a+/3)d:S 



^(p, a)[2, /3][r, a+^]c^S, &C. 



(6) 



-Lp-t^^-L ^^X^(XO 



Of these cases the first three require no special comment. The fourth case was for the 
first time solved by Professor J. C. Adams, in the Proceedings of the Eoyal Society, 
vol. xxvii., pp. 63-71, although the result had been independently found and published 
by Mr. Ferrers in his ' Treatise on Spherical Harmonics.' 

Professor Adams, as well as Mr. Ferrers, made the discovery of the value of (4) 
by an inductive process, and Mr. Todhunter showed in a su.bsequent number of the 
Proceedings how the proof could be thrown into a compact form. The present 
writer was, however, convinced that the integral could be found deductively according 
to the method described above, and the simple character of the result led him to hope 
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that the same method might be applied in cases of greater complexity. Accordingly 
cases (5) and (6) have been dealt with, but they have not been solved with the com- 
pleteness of case (4), the results being practically expressed in the form of series. 

The most interesting as well as the most important applications of the present 
method are connected with ellipsoids and eUipses. We have thereby obtained in 
series of harmonics the potentials of 

(1) An ellipsoidal shell. 

(2) A solid ellipsoid. 

(3) An elliptic plate of uniform density. 

(4) An electric current in an elliptic circuit. 

Of these cases (2) leads to the approximate determination of the forces expressing 
the mutual action between two solid ellipsoids or uniformly magnetized ellipsoids^ and 
(4) in certain cases to the forces between tw^o currents in elliptic circuits, the particular 
case of circular circuits being completely determinate. 

Similar results also hold for rectangular solids and circuits. 

The case (2) derives additional interest from the fact that it engaged the attention 
of Lagrange, who obtained the expansion as far as the first four terms. The same 
four terms have also been worked out in a very interesting paper on the potential of 
an ellipsoid at an external point, by Colonel A. R. Clarke (^ Philosophical Magazine/ 
1877, vol. ii., pp. 458-461). 

In regard to case (3), Professor Cayley, in the ' Proceedings of the Mathematical 
Society' for 1875, has obtained the solution in the form of an integral, from which he 
derives interesting properties of the potential depending on certain particular positions 
of the attracted point. The expansion in harmonic series would seem, however, to be 
practically more useful in determining the mutual forces between two electrical circuits. 

Throughout the following investigations, the method of treating Spherical Harmonics 
introduced by Thomson and Tait will be almost exclusively adopted. 



/x 



'""Fdl 



fh^ 



§ 2. Taking a point A on the axis of produced negatively at a distance a from the 
origin, let us consider the result of integrating 



C^n 



dS, or 

r 



z'^'udS 



over the surface of the sphere, where r is the distance of any point P from A, and u is 

the reciprocal of r. 

If we denote by the suffix 1 the fact of the operator being upon 0% and by the 

suffix 2 of its being upon % then the operator 
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dx dy dz^ 



upon the product, z'^ii, may be written, 



\ d? , d? , (d ^ d\n^ 

— +— «1_ I «1_ 1 V 

fZ^3 dy^ \dz^ dzj J 



that is, in virtue of u satisfying Laplace's equation, 

d^ ^ ^d dV 
_ I 2 ^ 

dzj^ dz^ dzr^l 
Referring now to the general result (3), and picking out the general term, we find 

47rR3^+2 i\ ^ 6^3i-« ^n 
___. O?^ . 

2i-\-l\ nl i—nl dz^ dz^ 
In order that this Diay lead to a value not zero, we must have 

or 

Substituting then this value of i, and differentiating z'^u, putting x, y, z=^0 after 
the differentiations, we get 

lb i ■ • 

2 
But u can be expanded in a series of zonal harmonics, viz. : it is 

Substituting this expansion in Hz'^'udB, and equating the coefficients of the different 
powers of the reciprocals of a to the values already found for them, we obtain finally, 
in the case where m and n are integers, 

2''ml — H— 1 
/x-P,(;?S=4:7rR^ : ^^^:^^ ,,..... (11) 




m + n-{' 11 



2 



It is obvious from the above proof that m-]-n must be an even number, and that n 
must not be greater than m. In all other cases in which m and n are integers, the 
integral must be zero. 
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It is obviouB the same method of solution as the foregoing will apply generally 
towards determining 



7(^) 



r 



d& 



where /"(z) is any function wliose differential coefficients are finite at the origin. In 
fact, we find 



ff/(z)P^c^S=Tr(2E)«+2. 

J J i 



n 



2n-\-ll 



'd' 



n 



E^ 



€^"+2 



+ 



E^ 



gn+4 



_dz ' 2(2n + S)*dz ' 2A(2n + 3)(2n + 5) dz 



+ &C. 



/(^)- (12) 



where, after the differentiations are performed, z is to be put equal to zero. 

A result practically equivalent to the above is given by Heine ('' Kugelfunctionen," 
second edition, p. 76), and applied to various cases where the functions can be ex- 
panded in simple series of zonal harmonics. The functions are 

Similar functions are also to be found expanded in Todhcjnter's "Functions of Laplace, 
Lame, and Bessel," §§ 48, 145-147. We shall now discuss more general cases of some 
of these functions. 



ff(l--/x^)-P,,rfS 



§ 3. The case now indicated will be discussed only for integral - values of m and n. 
Let us consider the integral 



c(x^-\-fy"' 



dS, or {x^+yYud8 



If we use the suffix 1 in any operator when it is upon (ir^+2/^)^ and 2 when it is 
upon II, we have as the general term operating, 

2i-{-ll\\dx-^ dxj \dy^ dyj dz^^ 
or, since 



d , d , d ^ d 



d d 
dx^ dx<^ d^^ d7]^ d^c^ drj^ 

d , d , ( d. d . d d 



^Vl %3 \^ll ^^1 ^?3 ^%j 



the general term above is 



2^• + 1 ! V \d^C^ dy\dr,;^ drij ^ dz, 
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that is, in virtue of ti satisfying Laplace's equation, 

47rE2^+3 ( d d d d d d V 
2^ + 11 " \d^j dT}-^ d^i dr]^ d^^ drjj 

The general term of the operator just found is 

2^ + 1! ul i—u\\d^-^drij \d^i drj^ d^^ drjj 

It is obvious that i—u must be even^ otherwise the differentiated expressions could 
not but vanish when x=0, 3/=0, z=0. We must thus have 

Substituting these values in the general operator, and retaining only the middle term 
of the last factor, which alone in that factor will produce a result differing from zero, 
we get, as the expression of the effective operator, 

4^j^3i+3 , m + nl Id d \^ [ d d Y 

4^ 



2m + 2^ + 11 m—nlnlnl \d^j drjj \d^^ di]^^ 



of which the last factor may be replaced by 

1 d^V 



TC Cf/^Q 



And now performing the operations indicated, we find as the result 

, . 47rE^^+^ , m + nl ml ml 2nl 1 

( lY 4^ — L___„„_ — __ — - _____ 



2m-\-2n-hll m—nlnlnl (X^'*+^ 

Since x^-^-y^^W' (1— /x^)=RV, and the term involving P^^ in the expansion of u in 
zonal harmonics is 



a^n+i-^ 3« 



we have, finally^ 



T/. o\ -r^ Id . T^o (— 1)**4»* m-^nlm\m\2%l , ^. 

(1— u^)^Po,,^S=47rR^T--^--^ — , , , ... * (13) 

J^ ^^ ^ 2m + 27^ + ll m—nlnlnl ^ ^ 



If (^ (v^), any function of v^, be expanded in the form 

Ao+Aiz/^+ . . . +A^+^z/2^^+^-^+ . 
then it may be shown that 

MDCCCLXXIX. 3 D 
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'1 



-^ ^^ 1.3.5 . . . (4»i+l)V »+ 4^ + 3 ^^'^^i + (4^ + 3) (471 +5)1.2 ^ ^»+i+'*«-y • ^^V 

In certain cases this series admits of summation. For example, let us examine under 
what circumstances JJ/x-^Pg^dS admits of being evaluated 



fji-''={l-i^)-i='Z 



, M+^-- (!+'>-' 



.j3n 



n\ 



Hence 



If we now write 



A --^ A 

A — \ / \^ / A Psrn 

[n ~\-l){n-\- A) 






^= n+l 

3 

y=2?^^-- 

we find for the integral J|ju,~'''P2»c^S, the expression 



where 



/# + 2)(Z; + 4) . . . (]c + 2n-2) 
, , «/3 ■ «(^ + lW^+l) , 

"^7"^ 7(7+1)1-2 ■^"••• 

r(7)r(7-/3-^) '"- 
r(7-ye)r(7-«) 

/ 1\ / l-/tf^ 

' ^+2) • • • ('"+^) 

^ Bertrand, ' Calcul. Integral,' 1870, pp. 495-6. 
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The value of tlie integral under discussion is accordingly 



(-1) 



\l -]€){?> -h) . . . {2n + l-'k) 



• • • 



. . . (15) 



and the above reasoning shows that h must be less than 1 : it may, however, have any 
negative value. 



(1-/>C^)^-*P2.C?S 



§ 4. To evaluate the integral now indicated, m and n being integers, we must 
modify our previous solution, for otherwise we should be led to infinite values of the 
differential coefficients at the origin. 

Let us consider the integral 



If V is a function which is symmetrical round the axis of z we may write this 

where c?5 is an element of arc of any section of the sphere of reference containing the 
axis of 0, and p^-{-z^-='K^, 

Now the expression (5) gives the form of solution in this case. If we put 



1 {'^ 
Y=: - (^+ip cos y\ff''d\\f 



we have, in fact. 



(l-/.2)«-iP2„(^S=; 



47rE2 



X 



jt+t; />2»(2+y/> cos t//)^# 



2^"+^«m + «.!m + «'! ]\dp dzj 

Expanding the quantities under the integral sign, and performing the differentia- 
tions and integrations, we obtain 



27rm^2m!29^! 



where 



g_^-^_K^^ + ¥) , <^--l)(^^ + i)(^- + ¥ + l) _ ^^ 



1.-2 



1.2 



1.2 



{n—m—^{n—m — \ — V) . . . ( — wi — J + l) 



1.2 . . . n 
3 D 2 
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First let n>m, then in S we have n—m positive and m negative factors, and in that 

case S becomes 

/ ^v 2n—2m\2m\ 

Next, let n<m, then S becomes 

, ,v 2m\m—n\ 

Substituting these values of S in the expression for the integral v^e find 
(1) if^^>m 

{-iy'2ir^W2m\ 2n\ 2n\ 2n--2m\ 



23»H-4?i^^ _j_ ^^ I rifi_qji\ ^l ryy^l ^j ^ 



(16) 



(2) Jin<m 



(-l)^27rm^2m! 2m! 2n\ m-n\ 
2^'''+^''m-{-nl2m—2nlmlmlnlnr •*••'•'■• V / 



ffQ.fdS.ffcgd S 

§ 5. The results for the two integrals now stated are, if a be the colatitude of the 

pole of Q^-, 

47rE^P,-(cos a) .......... (18) 

27TU\j-i){j+i+l)Pi{co8 0L) . ....... (19) 

where, in the second integral, y—^ must be a positive even integer. 



G-ENERAL THEOREMS IN DIFFEREKTIATIOlSr. 

§ 6. Before we proceed with the remaining cases, it will be convenient to state and 
prove various theorems in difierentiation."" 

Theorem i. The general operator (7) upon the reciprocal of r may be made to 
operate, instead, upon a homogeneous function in x, y, z of the i^^ degree ; this 
effect being expressed by the relation 

d} 1 1 d} 

^ ^ dh^dhc^ . . , dhi r'^r'^^ d\dh^ . . . dhS W?; ««••«« \ ; 
where the pole of Q^ is in the direction of r. 

* Theorems i. and ii. were given in a paper by tlie author in the ^ Messenger of Mathematics,' No, 73, 
1877. On account of the brevity of the proofs, and in order to secure completeness, it is thought best to 
reproduce them here. 
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Theorem ii. relates to the differentiation of the reciprocal of r^*^"^^ and may be 
stated thus : 

(_l)i-. <tl L_=2^! *Z 1^-/1^^ . . . (21) 

^ d\dlfh^ . . . d'hi^^T^'^^'^ 2cr! d\dl\. . . dhi_\ dyiF j ' ' ' ^ ' 

Theorem iii. 1£ f (y, --, -y) be any homogeneous operator of the #^* degree 
operating on a homogeneous expression <^ (x^ y, z) of the i^^ degree, then 

J. ( d,- d d\ , , . ,fd d d\ ^ , . . . 

■^U'J/^J^(^'^'^) = n5^'^/&j^^^'2^'^) (22) 

This theorem is obviously true, and though stated for only three quantities, x, y, z, 
is true for any number. 

Theorems i. and ii. are proved as follows: Let OP=r, QP=p and oPQ=it—0^ 
where Q is any point near to P. Then 



OQ r r^ 



fl+ . . . +(-l).JjQ,+4o. 



Keeping r fixed let us perform the operation (7) on the two sides of this equation 
and then put p=0, in which case OQ becomes r. It is to be observed that with a 
homogeneous operator like (7), it is immaterial what point is taken for origin of 
coordinates. If then we take P for origin when we are dealing with the right hand 
side of the above expansion, we see that the first i terms will disappear by repeated 
differentiations, and the terms beyond the i-^-V^ disappear in consequence of the 
zero value of p. The theorem stated therefore follows, the substitution of r for p 
making no difference in a result from which x, y, z are finally made to disappear. 

Theorem ii. may be proved in the same way, if instead of the expansion for tlie 
reciprocal of OQ, as above, we take the expansion for the reciprocal OQ^"^"^^, found by 
differentiating the above series a times with regard to jlc. 



DIFFERENT FORMS OF TESSERAL A¥D SECTORIAL HARMONICS. 

§ 7. In several of the following investigations frequent use will be made of different 
forms of the zonal and tesseral harmonics. With the view of classifying these various 
forms, and of bringing the various expressions for the tesseral and sectorial system 
into harmony with the corresponding expressions for the zonal, a proof is here given 
which will be a direct illustration of the foregoing theorems. 

The most important forms of the zonal harmonic P^ are 

1 d^ 

2H^' da ^ (A; 
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Now by definition 






e e 9 « 6 o 



{lJi-{-jv COS xjjydijj, - 

'^'^ 



TTJ 



IT 



d^ 



(/J. -hjp cos 'fy 



+1 



a e c G c 



*■(*■-!). ,,-_o, 3 



1^^+ • 






-^ ! dz\d^ diri)\/z^ + ^?7 



• (B) 



(C) 



9e»» 09 0$ 



(D) 



By Theorem ii. this is 



-(-ly 



+0-^ 



2r' 



'i(T 



2^H\ a\ 



rp 



^•+1(^^ + 7?^) 



dr-"" 1 



2H\ 



-2 cos (T(^v''-j r 



d}-""! . d^Q 

dfJb^ 



dz 



We observe that, in this result, the axis of the harmonic Qi is the hne OP, as explained 
in the proof of Theorem i. We must therefore change the axes of coordinates, 
making the new axis of z coincide with OP. 



d-V, 



If r''"'^-^', when expanded, becomes 
we have to consider the effect of operating upon thi^§ with 



d . . d\"^ 
cos a— — sm a-~ 
dz dp) 



There will result 



that is 



^ — cr! (A cos^ ""a— Bcos^ "^ ^asin^ot+ . . .) 



^ — cr ! 



\ -^ 



djjF 



where \l is now the cosine of the angle between OP and the axis of %, nence we 
obtain the first and second forms of the harmonic corresponding to the forms (A) and (B) 
of the zonal harmonic, viz. : — 



(i'^)='^'2coso-.^v'^-7;;. 



2^ 



dpf- 



a s e 



99A&s>6 oe« 






2'+'^il i 



^— (T, 



2(2^-1) ^ -r • ■ -j ■ 



9 « 



• («) 



• m 



Let us now consider the definite integral 
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,Tr 



dcf) 



TT 



p cos^ ^ + S' si^^^ ^ 2 VM 



• e«»«»»» x^^y 



This result will still be true if for p, q we substitute the unreal quantities c+re^^', 
c+re"^^', as may indeed be shown by direct integration. 

Let us diiferentiate both sides of the equation (a), or times with regard to p and a 
times with regard to g'. We find 



2cr! 



•1 

2 



sin^'^ ^ cos^"^ <p 



(p cos^ S^ + 2' si^^ c/))^"^"^^ 



C^^: 



TT 

'2 



L3 . . . (2(r-l) 

2c. J (^^)cr+| 



Making the above substitutions for p, q we get 



2^al 



2f} 



sin^'^ 2(f) 



1.3 . . . 2o--l ttJo {G-{-r(e^^ cos^ (j^ + e'^J sin^ ^)} 



2(r+l 



C?(j6 = 



2o-+l 



^c^4-2cr//,4-r^^ 



1.3 .. . (2(7-1) c^^+iV^Z^- +...+( 1)' 1 



'^\i-o- ^0-p 



d/ju^ 



I • » • 



the last line being derived from its predecessor by a differentiations with regard to /x. 
And now expanding the expression under the integral sign in powers of r, and 
equating the coefficients of r^""% we find 



d'^Ti i + crl al 



2^- 



d/jL"^ 2a\ i—al tt. 



2 




{e^-^ cos^ (j^+e"^-^' sin^ ^)^-'^ sin^^ 2^d<f> 



On putting 2(f)=z\jjy and substituting the expression just found in equation (a), we 
finally obtain for (i, cr) one of the forms corresponding to C, viz. : — 



i-i-al al If" 

-y-pj-2 COS cr^i^^ - (/>t+jV COS i//)'""^sin^^ i//(ii// . . 



. . . (c) 



There is obviously another expression, corresponding to the second form of C, which 
is easily written. It is to be found by expanding, as above, in powers of the reciprocal 
of r. 

The expansion of the binomial inside of the integral (c) and subsequent integration 
with regard to i// give the form corresponding to D, viz. : — 



2^Hl a I 



2 cos cr(j)v''\ fj} ""' 



(i-o){i-€r-l) ._^ 



4(a-+l) 



F 



•2, .2 



V^+ . . 



• . 



. . {d) 



Integration of the Product of Tim Harmonics and Proo/ of Laplace's Coefficient 

§ 8. Let IIS take two points A and B at distances a and h from the origin, A being 
on the axis of z and B at an angular distance a from it. 
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Let P be any point on the sphere of reference, and let ns consider tire integral 

AKBP 
By (3) we may write as the result of this 



2i ~h 1 \\dx dy dz ] AO.BO 

If we use the suffix 1 to denote differentiation upon AP only and 2 upon BP only, 
the above becomes by virtue of Laplace's equation 

2t + 1 ! \dx^ dx^ dyi dy^ dz-^ dzj AO.BO 
or, what is the same thing, 

^^^ o2i + Tl \"^ drj.'^cl^., dTj^^dz^ dzj AO.BO 

Now on consideration of Theorem i, we will see that in expanding the operator just 
found, the only terms which will lead to results not zero will be those which, so far as 
the operator with suffix 1 is concerned, are of the form 

^m gm fp 

d^i d7]i dzi 

and this by virtue of Laplace's operator 

cV' ^ d d ^ 
— j-4 =iO 

dz^ d^i drji 

d 

is expressible in terms of — only. Hence we conclude that the only effective terms 
may be found directly as those not containing a power of v in the expansion of 

4:iTW%7:T-rT,{ — r'7 -, -z — 1-2- 



^2^4- 111 dz-^ v^ dz^ dz^ dz J AO, BO 



» • 



CO 



f/tty» a 



4irR^S(— iVrT-T-Tt '^y~"~~T^ Tnwi 
0^ ^2t.+ ll\ dz^ '^ dz^j AU.Bu 






xxexiCio 



'2^ + 11 i\ i\ dz^ dz^ AO.BO 



d^ . -ri9v< "-P^^* Pi (cos a) 



JAP.BP"~" "2i+l c^^+W+i 
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If now we expand the reciprocals of AP and BP in harmonics we arrive at the 
well-known theorems 



^hiQjdS=0 



•3 



4:7rW 
P^Q.(iS=-^— Pi (cos a) (23) 

The last result, when properly considered, can be made to yield the integral of any 
two harmonics. For it can be easily seen from the foregoing work that if p^, p^ are 
any two vectors drawn from the origin, then we have this very remarkable formula 
for Qi 

2i\ Qi=p/+W 2^ -7- 7^+-^- -1 — H-T" T' "^ 

' -^ ' ^ \dx^ dx^ dy^ dy^ dz^ dz^j pip^ 

On account of the operator being an invariant we may suppose the axes to be any 
rectangular axes whatever. 

By what has just been shown we may throw this into the form 

/ d ^ d \3« 1 

2ii Q.=p/+V2'''H-im3^--#-) -^ 

^^ HI H2 \ ;y ^^^ ^ a^j p^p^^ 

where v^ stands for 

d d 

» 

d7)^ ' dif}^ 

The general term of the value of 2i\ Q^ just found is 

i — a\ i-{-(T\ v''^'' dz^ dz^ p^ p^ 

6.6., 

p/+^P3^+X — 1)^2^1 d<' d-"" d'-^' 3^ 1 
i—a\i-\-a\ drj^ drjc^ dz^ dz^ pip^ 
or 

i—crl i-{-a\ drJY d^^ dz^ dz^ p^p^ 

It is obvious a term similar to this can be obtained with only ^^ put for ^^ and 17^ for 
7)1^ by merely changing the sign of or in the first of the last three expressions. Adding 
the two terms and recollecting the definitions of the harmonics we may write the result 

^^^^^^^ 22--i{(io-)(icry+ft<T]y, (xj} . .... . (24) 



i — (T\ 'i + cr! 



Hence Laplace's Coefficient as given in § 1. 

Reverting now to equation (23), let us take two fixed poles, and let Q-, Q/ be 
expanded according to Laplace's Formula; we shall then be at once led to the well- 
known surface integrals of two tesseral harmonics. 

MDCCCLXXIX 3 E 
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The discussion entei^ed into in this article is necessary for the treatment of the 
product of three harmonics to which we now proceed. 
It may be remarked that the expansion of the operator 

el d d d d dV 1 

dx^ dx^ dy^ dy^ dz-^ dzj p^p^ 

is, according to the above investigation, 



i\i\ dz^dz^ * * i-\-(r\ i—al \d^i dr}<^ d^^ dr}J\dz^ dzJ ' ' J p^p, 

(See also Thomson and Tait, p. 157). 



General Method of Dealing loith the Integral Qj^Q^Qr . . . dB, 



J .) 



§ 9. In discussing the general method of evaluating this integral it will be con- 
venient to confine ourselves to the case of three harmonics, though the first steps of 
the reasoning will apply to any number. 

Let A, B, C be any three points whose distances from the origin are a, b, c, and let 
P be any point on the sphere of reference. Then by § 1 

. T^n^ R^^ fcP , d^ , d^V 1 



AP.BP.CP 2i-^ll\dx ' dy ' «^^/AO.BO.0O 



2^+11 ^ "^ ' ^' ' -^^'^ AO.BO.CO 
where %3 stands for 

d d d d d d . 

CvXn CoXo Cv U q Cv U o (l/Zn CvZo 

the sufiixes having the same signification as in the previous article. 
The general term of 2i-\- 2 dimensions in B may therefore be written 

2^+11 Vl /.! z;! ^ ^s^'si'^^s^o BO.CO 

Now the general expansion of each of the operators in this expression was found in 
the previous article. 

If therefore we adopt a somewhat more convenient notation we may expand %3^ 
%i^> ^'i/ 11^ ^ series of terms of which the type is 
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o,;a. J.S ' 2X12luI 2v\ 

X + l\ X—ll fM-hml fjL—ml v + nl v—n\ 

I d} d} ^ d} d^\/d"' d^' , d'"" d'^'Xf d'' d"" . d'' d"" 
■v ( I ~— ]( -J 11 -4- — — - — 

\^?3 ^% ^^fs dnj\dh dVi d^i dvj\d^i d% d^^ drj^^ 

^+v—}n~n ^v\-\—n—l ^K+fx-~l~m 

dzj dz^ dzo^ 



where the values of I, m, n range between and X, /x, v. 

The general operator just found we shall denote by the symbol tt {l^ m, n). 

§ 10. If we expand each of the reciprocals of AP, BP, CP in harmonics we shall 

arrive at the general term of 2i-\- 2 dimensions in R in a different form, viz. : it will be 



R3^ 



^i»+i^^+i^f+i 



^p^q^r^^ 



In order therefore that the general term as found by the work of the last article 
may correspond with this we must have 

l^"~|~* A— — (j 

The quantities X, /x, v, as depending upon p, g, r, are therefore perfectly determinate, 
and the equation expressing the identity of the general terms may by Theorem i. be 
written 

%%QrdB=^^^2;^^^;~^^tt^^^^ . . , . (25) 

On examining the equations determining X, /x, v, we see that in order to their being 
positive, P+9'+^ must be an even number, and no one of the three quantities p, q, r, 
must be greater than the sum of the other two. 

If these conditions are not satisfied, then JJQ^Q^Q;.c^S=0 for integral values of p, q^ r, 
wherever the poles of the harmonics may be. 



JL ^JL rtJL j'Oi/kJ 



§ 11. Let us now suppose that the three points A, B, C are in the axis of z. Then 
the harmonic Q^ becomes P^, which we will suppose expressed by the series D, In like 
manner P^ and P^ may be similarly expressed. It follows that, in selecting from the 

o IS >4 
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general operator what terms are practically effective, and discarding those which are 
not, we need retain only sucli terms as are, for example, of the form 

d^ d^ dP-^"^ 
d^^ dr]^ dzj 

Now the terms of this form are obviously got by putting lz=zm^=-n, =cr say. In that 
case, when we omit inoperative terms, tt becomes 



>6cr-- ^'+1. 



2\\2fji\2v 



X 



X + cr! \ — G-\ jjb-^-al /jb — al v-\-a\ v — al 
d d d d d d Y dP-^^ d^-^^ d'~^^ 



^d^^ dr}^ d^^ dnr]^ d^^ dri^j dz^ dz^ dz^ 



that is, in virtue of Laplace's equation being satisfied, 

, 2 2\\2tJi\2v\ dP d^ d>' 

l—\Y 

^ ^2^ \ + cr! \ — cr! /^ + cr! fju—crl z^ + cr! v — al dz^ dz^ dz^ 

Taking all possible terms of tt, such as that just found, and performing the operation 
directed by them, we have for jjP^P^P,(iS, the expression 

where cr ranges from to the smallest of the quantities X, /x, v. 

Now the series within brackets is clearly the coefficient of x^^, y^^, %^\ or of 'xr-'^% 

or, again, the same series is also obviously equal to 

plqlrl ^ 
2\\2iJi\2v\ 

where K is the coefficient of x^"-' f"^ z^''^ in 

{—lY{;y-'Z)^\z—xY^{x—yY'' 



tliat is, in 



(_i).<^-r*v-( ^y- ^/tf{ ^\- ^m V=- V; 



yl\^x ^ zj \ ^ y ^ x^ 



%> 



If we now write herein 

zzTzxe'^^'^ye'^^ 
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we see tbat K is the term not involving the cosine of an angle of the value md-^n^ in 
the expansion in such cosines of W, where 

W= (— 1)^' 2^ sin^^^ sin^'^c^ sin^"(^— ^) 
so xnatj 

Now expanding sin^''(^— ^) in powers of sines and cosines by the binomial theorem, 
we observe that the odd terms in the expansion disappear on integration, so that the 
value of K is given by 

^h ^{\- tf^^^^!^^-~ sin^^+^'d GOB^^-^^d sin2^+2v~2^^ ^^^2.^ ^g ^a 

On integration and reduction this becomes 

/ i\,;n-~2t. 2^' ^ 2\-]-2%\2iJb + 2v — 2n\ 

( iVo-~2 . ^^'' 2\\2iM + %v\ 1^ 
X + vl fi + v\ X! z^ + z^l v\ 

q^-l, (2X+1> (2X + 1)(2X + 3>(^.^1) .n 

"^ 2/^ + 21^-1"^ (2//, + 2z;--l){2/^ + 2z^ --3)12"^ 

This hypergeometric series is capable of summation by a method similar to that 
given by Bertrand, ' Calcul. Integral/ 1870, pp. 495-496. If we put 



that is, 



wh 



ere 



S becomes 






1 _L^^_L ^(^"^^) K^""-^) _i k 

'r(a)r^+i)Jo(i+2/)^+''+' 

r(«+^+i) T(fi+i) 
T(^+i) r(x+/i+i) 

X + /U.' /x! 2^ + 211! 
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.On substituting this value for S in tlie value of K we find 

Finally^ performing tlie operations witli regard to tlie z^s, we find 

47rR^ il i\ 2X1 2/^! 2pI ^ . 

^P^s^rdb=^-^^-^^^j^^^^ (ib) 

This agrees with, the result obtained by Professor Adams. 

§ 12. Before leaving this case we may remark that, so far as the foregoing proof is 
concerned, the poles of all the three harmonics are not necessarily coincident. The 
proof will bold when two are coincident and the third is at an angular distance, say a, 
from them. In those circumstances the value of |JP^P^Q;.(iS will be the result found 
above multiplied by P,.(cos a). The same is, of course, true of the integral ||Q^Q^P,.(iS^ 
provided the two Q's have the same poles. Now let Q^, Q^ be each expanded accord- 
ing to Laplace's Formula. Then we will have the curious result that any zonal 
harmonic of degree r can be expanded in products of harmonics, zonal and tesseral, of 
degrees p and q, the quantities p, q, r being, of course, subject to the foregoing 
restrictions. If q is less than p the number of terms in such expansion will be 2g+ 1 



(p, a){q, fi){r, a+fi)dS 



§ 13* Reverting to the general expression for the operator ir (I, m, n) given in § 8^ 

let us leave out for an instant the numerical multipliers and the operators with regard 

to the z^8 and let us multiply out those with regard to ^ and rj. Let us suppose Z, m, n 

in descending order of magnitude, and for the sake of brevity^ instead of writing the 

■J 

differential operators, let us write only the characteristic letter, e.g,^ instead of yr- 
write ^j. We find 

d d 1 d^ 

Now since any product such as j^ j-~ can be replaced by — . _ ^ wa see that each of 
the lines just written corresponds to the sum of tesseral harmonics^ in such manner 
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that if we suppose the differentiations performed by the whole of the operator, and if 
one of the resulting terms is 

w^here the suffixes indicate that the angular coordinates of A, B, are substituted 
in the harmonics to which they are respectively attached, then we have in all cases one 
of the quantities a, /3, y equal to the sum of the other two. This we see from the 
expanded form of the operator just given, and it is a result which was to be expected, 
for supposing Q^, Q^, Q,. severally expanded according to Laplace's Formula we should 
have in the surface integral a series of terms of which the type involves the integral 



r2Tr 



COS <x^ cos fi(j> COS y(f)d(f) 

and this integral vanishes unless the above conditions be complied with, the case where 
a4-iS+7=0 being of course out of the reckoning. 

In accordance with the explanation just given we may now put 

m — n=a m'^7i=a 1 

or I 
/+9^=^ l^7i'=/3 y (27) 

and /. l-}-m=a-\-/3 J 

If then we substitute in tt {I, m, n) for I and m the values ^+n and a— ti, we see that 
n may range in numerical value between and the least of the integers X— ^, /x— a, v, 
and may be positive or negative. The operator tt, in fact, becomes 



/ 2 \«22('^+P)-~'' ^^ ' ^^ • ^^ • 



X + ^ + ^l \—j3—n\ fi-\-a—nl fi—a + nl p-\-nl v—nl 



K^^i ^^3 ^Vs ^^Vi ^V2 ^^fs/ ^^% ^% ^^. 



3 



We have therefore now to consider a series in which any term is the product of (1) 
an invariable power of 2, (2) an invariable operator, (3) a variable coefficient whose 
parameter is n. If we denote the first two of these factors by the symbol ot, then the 
sum of all such terms is t»yK, where 

and W is the coefficient of x'^^y^^z^'' in 
that is, in 



400 MR. W. D. NIYEN ON CERTAIF DEFINITE INTEGRALS 

The expression upon which the operator m has to operate may be found by expressing 
Qoj &c., by Laplace's Formula, and using the series D: it is as follows : — 



a\ 'p—al 

If we omit irrelevant terms this becomes 



ol\ /3l ci-\-/3l p — oil q—/3l r — a— y8 

-{P> ^)fe ^][^^ ^+/31~[l>. ^Ife ^)^ [n ^+/3J~[b ^Ife ^Mr, a+/3).} 
The result of the operation is accordingly 

Turning now to the integral HQpQqQrdB, let us expand each of the Q's according to 
Laplace's Formula. We thus find for the general term 

2'^(^+^^-yi p\ q\ gl rl rl 



p-\-al p—ul q-^fil q—/3l r + a + /3! r—a—^l 



r 



X 



J J 



{(p> '^)^ip> ^)+[p^ ^l[p^ ^]}{(^^ fi)h{q, 0)+gj ^lfe^]}{(n «+/3).(n ^+^) 

Now, by the formula of § 10, when we come to equate the two results just found, we 
shall get 

2^i'^-^^-Y. ,1 rl _ CC , ^ i^ il 2X! 2f.l 2.! 

or finally 

where W is defined above (28) as the coefilcient of a term in a certain product, or 
as a series whose terms depend upon a single parameter. 
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If we denote by U the value of the expression on the right hand of the equation 
just foundj we have further 



JL 41 -L f/L f'-i- g\A/hJ 



§ 14. We may remark in regard to the case which we have now to examine, that 
since P^P^ ma.y by the foregoing results be expanded in a series of the form 

where the value of A,^ determined by (26) is 

2r— 25 + 4m + l r-j-ml r-\-m\ 2ml 2r — 25 + 2ml 25— 2m! 



2r + 277^+l 2r + 2m! ml ml r—s + ml r—s + ml s—mA s—ml 
the above integral is reducible to the computation of a series of s terms of the form 

r t 

The last integral is to be found by writing in (26) 

2X=-— p-f-g4"^"--"^+2w^ 

2j^=p+g— r+5 — 2m 

We are thus led to a complete though practically laborious evaluation of J|P^P^P^P^cZS. 
It will be observed that the result is 2jero unless p-^-q-^-r-^s be an even number. 

With the view of finding* out how far the method of this paper is applicable in this 
case, and what difficulties stand in the way of its general application, we will now 
briefly apply it. It will be convenient to expand the operator in a somewhat different 
manner to that pursued in the case of the product of three harmonics. The operator 
was then expanded in a form which would render it useful for application in the case 
of the tesseral and sectorial system. If, however, the poles of the harmonic are all in 
the axis of z^ a much simpler mode of expansion may be adopted. 

When the product of four harmonics is under consideration we have to discuss an 
operator of the form. 

MDCCCLXXIX. 3 F 
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d d d d\^ 



, xJbJb-x CtJf^g Ojth o 



s XT "T" XT "T ^T" I + 



d d d 



dV- 



9, 1 



dxj \dy^ dy^ dy^ dyj \dz^ d%^ dz^ dzj _ 



or, what is tlie same thing, 



±+±+l,+A\\,, 



d 



dz. dzcy dzo d% 



d d d 



■4/ 



d d d d 
,d^i ' d^^ ' ^^3 • d^J\dr]^ drj^ drj^ drj^ 



)} 



If this operate upon the reciprocal of pip2Psp4,^ where these distances are measured 
from four points A, B, C, D on the axis of ^, we have four relations of the form 



d^ . ^ # 

— •4-4 

dz d^dT] 







Further, on picking out the general term^ we observe that, to be effectually 
operative, it must consist of products of the form. 

d^ d''' d''' 
dz d^ dr) 

and this form of operator may, by the foregoing relation, be expressed in terms of 
d ^ 

It is obvious, however, that the sum of the terms so modified can be found directly 
from the part of the expansion not containing h'^ of 



d d d Cl\^ /7 . 7 , 7 . 7 x/1 ^^^ . 1 ^^^ , 1^^ , 1 ^^ 

^_^^~— -L-p-|_-- «-./^ J_^ _j«^ _j_/. \ J^ _j_ ---[-- ~- 

dz^ dz^ dz^ dz^l n j- • ^ ^ -t/ y^^ ^^^^ j^ ^^ /^^ (^^^ j^^ (^i^^^ 
or, what is the same thing, from the part of the expansion not containing ^'s of 



d d . d ^ d\^ f T d . T d ^ ^ d . T d 



S^Z-t CuZa CvZ, 



^%. 



*^"*""' ' ' '^" iiii 111 III 1 1 — *~— ■ ■ ''^ ■«<■•■■■■ ■" 



CLZ-t CvZn Ct'/vq Ka/ZaI \i\i-\ CvZ-i fVcf Cv/i-g /i/Q LvZo '^d, ^^Za 



1 ^'^^l 



ri i^ 



.3 ct.^^3 ,v^ vv^^^ 



that is, of 



<-^)'{(Vl-V: 



-^ ) v~ 7"+5 similar terms 

hj dz^ dz^^ 



If we expand this the general term will be 



(-1)^ 



i\ 



X! [julvl II ml n\ 



X(^,"»-^s)^^fe-^l)^'^(^l~^.)^1^4-^l)^^(^4-^.)'^1^4-^3)^'^ 



X 



1 d \P/ 1 d Vifl d V/l d 



ki dzjj \k^ dzJ \/^g dzJ y% dz^ 



where 






=p 



"^ 



>' 



(30) 
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The quantities X, fi, v, I, m, n are therefore indeterminate, and there will be a series 
of terms having the required indices. 

The coefficient of the operator given above ma,y be put into another form. Let u.s 
consider the expansion of 

first in powers of a, b, c, . . . and then in powers of k^, k^, k^, , . . We then see that 
the coefficient of hf ^/ k{ h^ in (^3—^3)^^ . . . (^%— ^3)'^'' is equal to 

2\\ 2tjL\ 2vl2ll 2ml 2^' 

r: y^ 

pi ql rl si 

where W is the coefficient of a^^ ¥^ c^"" (P^ ^zmj-^t ^^^ 

We have, finally, 

T> "o T> T> 7Q— / ^K-^'n-UHl 2X1 2fil 2v\ 211 2ml 2^^! 

J:^F^r,r,db— (- Ir^T-^ S XI fj^l vl II mAnl 

where t denotes the sum of all such values of W multiplied by their respective 
coefficients corresponding to values of X, jju, v, I, m, n determined by equations (30). 

This investigation serves to exhibit the peculiar practical difficulties which beset the 
problem of integrating products of harmonics over the sphere of reference, if more than 
three harmonics be considered. 

POTENTIALS OF ELLIPSOIDS. 

§ ]5. Let p be the perpendicular upon the tangent plane at any point of an 
ellipsoid, and let us consider the integral 

» 

taken over the surface. 

By the theory of corresponding points this integral may clearly be thrown into the 
form of an integral taken over the surface of a sphere of radius R, viz. : it is 

the value of which, by what we have already shown in (3) is 

27raoc- — -7— 

V 

where 

3 F 2 
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Now let V be any function of x, y, z, which can be expanded by Taylor's Theorem ; 
then 



J J 



taken over the surface of the ellipsoid is 

2iTabc — — — V 

V 

or^ in series, 

w3 r-rSi 



where 



o I Zi% -f- 1 ! 

clx dtfi dz ^ 



and after the differentiations are performed, x, y, % are to be put equal to zero. 

§ 16. If V satisfies Laplace's equation, then the operator v^ will remain unaltered 
if for c^, W, 0^ we substitute a^+^^ 6^+^. c^+e. It follows, therefore, that the average 
value of V over any ellipsoid as measured by 

• .0 • > • • • • • > \%J> X, t 



o 

O 



volume of ellipsoid 



is the same for that ellipsoid as for any ellipsoid confocal with it. 

This theorem/"'' which may be regarded as the anologue for ellipsoids of the corres- 
ponding theorem given by Gauss for spheres, is due to Professor Charles Niven, who 
also showed that, if V be due to attracting matter E inside the surface of the ellipsoid, 
then the expression (31) becomes 

Ep df 

1 
§ 17. As a particular case, let V=" where Q is any point outside of the ellipsoid. 

J, AltJix 

JjQP 

is the potential at Q due to matter of density p coated over the ellipsoid. The 
quantity of matter is ^itabc and the potential due to it at Q, according to a well 

known result is 

d-^ 



^00 



•'e 



-v/((^2 + t)(6^+i^)(c2 + ^) 



or 

where a is the ellipsoidal coordinate of Q. 

* Math-ematical Tripos Solutions for 1878. 
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If P' be any point inside the ellipsoid we then haye ultimately when P' coincides 

with the origin 

^v_g-v 1 
a= 



QF 



Bnt, if y^ g, h be the coordinates of Q, 



Hence 

d 1 d 1 



dx QP' df QP' 

dx qv'"' df qv 



&c. 



We may therefore take V ^ to stand for 

"^ df'^^ dg'^'' dli 

and in that case we may put x, y, z zero before differentiating. If then OQ be 

denoted by r^ we have 

aV-g-v 1 

06= .......... (32) 

y r ^ ^ 

or, in series, 

i^+b'+ ■ ■ ■ +dr, ^"+ ■ ■ ■)l (33) 

If we multiply a by 27TpabcBd we have the potential of an infinitely thin shell 
bounded by two similar and similarly situated ellipsoids. When a shell is hereafter 
referred to, this kind of shell will be meant, unless otherwise specified. 

§ 18. Since the whole ellipsoid may be divided up into shells, the potential of the 
ellipsoid at Q is 

P^qp^ (34) 

where p is the density at any point in the shell whose semi-axes are a9, bd, c6. It 

will be observed that, since pSd=OSp, - is an element of volume, so that (34) may 

be written 



f 



^ pdxdyd^ 



If we take the series (33) we find for the potential of the ellipsoid at Q the series 

'\7rpabc0^(l-ir^^ V^^- . . . +-^ V^^+kc^^dd .... (35) 

I \ 6\ Jt + ll It ^ ^ 
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Tliis can be integrated if /> be a function of 0, 

In the particular case where p is constant and the mass is denoted by M, we find 

-1. tJL I -1. I •^ 4-« J y I • • • 1 • I /-i # iTk r^ t , -t t I •••I • • • * • \ O \j I 

\ * 5 3! ' ' 2^ + 3 2^ + 11 Jr ^ ^ 

§ 19. The result jnst formed exhibits in a very simple manner Maclaurin's 
Theorem, that the attractions of confocal ellipsoids at external points are proportional 
to their masses. For, as has been already pointed out, the operator V^ is unaltered 

/(P (P (P\ 
by the addition of ^ [tA-^ — ^""^ )• 

The series (35) also shows in what direction Maclaurin's Theorem may be 
generalised. 

(i.) Let pz=zf{9). We then see that if there be two ellipsoids which would be 
confocal if they were coaxal, and if the matter in them were arranged in layers 
similar to their bounding surfaces according to any specified law depending on the 
parameter of the layer, then the attraction of such ellipsoids at external points whose 
coordinates referred to the axes of the ellipsoids are equal, are proportional to the 
masses of the ellipsoids. 

(ii.) If we multiply the value of a given in (33) by 



27rv'(a^-<^)(6'^--^)(c^-</>)F((^)d<^ 
and integrate between the values (jy^ and ^3, we find on the right hand side 

An\'y{a?-4>W-<i.){c^-4>m<f)d<i>x{l-¥^+ . . • +2^,+ .■.)l . (37) 

in which the operator is independent of <^. 
On the left hand side we have 



{^\rr^{a^- ^) {W- <^) (#- <^) oF {<f>) # 






4> 



J J 



yi 



where ^ is the perpendicular from the centre on the tangent plane at the point x^ 3/, % 
to the confocal, passing through the point of parameter <^. 

Hence, if the ma-tter between two confocal ellipsoids be affected with a density 
varying as F((^)p^, the equipotential surfaces will be confocal ellipsoids, as may be seen 
from the expression (37), which also gives the ratios of the attractions due to different 
sneiiSa 

(iii.) The theorem is true for any value of V satisfying Laplace's equation. 
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§ 20. If we put a=::=& tlien V^ becomes — (a^— c^)--, and the potential of a solid 
ellipsoid of revolution of uniform density becomes 

3-5 T^ '^b-iX T^ 1 7-9 V r^ /"+"•••; * * ' ' ^"^^^ 

§ 21. The expansion involving surface spherical harmonics of the potential of a solid 
ellipsoid of uniform density is easily derived from (36). 
In virtue of Laplace's equation, the operator V ^ becomes 

which we will abbreviate to 

clx dz 



Now by Theorem i. 

dx ^ dz J r ^ ^ T^'+^\^ dz -^ dxj 



Since there is here no operator in regard to y we may put ^=0 in the expression 
for Qgi'?'^^ before differentiating. In that case ^=^7)=^x, and the expansion will take a 
comparatively simple form. 

We take, according to custom, the axis of z for the axes of the zonal and tesseral- 
sectorial system. The expansion of Qg^r^* then becomes, by Laplace's Formula 

04a— 1 9V! 9o'\ 



4 






+2(7! 2^-2(7^'^ '^'^>' V^ "" 4(2(7 + 1) ^ "" +...)+ &c. 

The general term of the expansion of the potential of an ellipsoid, may therefore be 
expressed thus :— ™ 



O 1 

O J 



where 



, (i—or)(i—(T — l) 2<r+3 1 o- o ^ />-> . ^ n \ 
+ ^—42 -2o^2^9''-'-r^''-^o.). . . . (39) 
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§ 22. It may be observed that the symbohcal form of the operator in the case of the 
potential of a shell leads to a symbolical form for the potential of the ellipsoid, viz. :— 

3M v(^v^+f "^) -:^^ 1 

Now let there be a second solid ellipsoid of mass M^ and of semi-axes a^, ¥, c^ and let 
it be placed in a position where the coordinates of its centre referred to the axes of the 
first ellipsoid are /, g, h, and its axes are inclined at direction cosines (Ij m^ n^, 
(4 m^ %), (?g ?% %). Then if we put 



wi: 



here 



■^ dm at/ dz 

d -1 d ^ d . d 
dx ^df ^clg ^dh 

Cv 7 I Cv Cv 

'^^^/^ '^ '% ''^' 

d J d d d 

dz^''^dr^^^^g^%E 

a double application of the reasoning of §§ 1 6-1 8 leads ns to the conclusion that the 
exhaustion of the potential energy of the two ellipsoids due to mutual action is 

9^^^^^V(e ^-^e^'^)—e^ + e~-'^ Vx (e^^ + e'"'^0— e'^^ + e"^^ 1 .... (40) 

where r is the distance between their centres. 

If we turn the second ellipsoid through an angle B6 about the line (X, /x, p) passing 
through its centre, we have 

Bmj^=:{Xni — ph)^0 r ..»...,.. (41) 

with similar increments for the other direction cosines. Hence the expansion of (40) 
in harmonics leads to an approximate determination of the forces and couples repre- 
senting the mutual action between the two ellipsoids. 

It is obvious that a similar investigation will apply to the determination of the forces 
and couples between two magnets in the form of ellipsoids uniformly magnetised. 

It may be interesting to notice that the foregoing method of expansion (36) shows 
that for points at a considerable distance the potential due to a solid ellipsoid is the same 
as if its mass were distributed as follows :— Two-fifths at the centre and one-tenth at 
the extremities of each of the axes« 
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General Remarhs on Ellipsoidal Surface and Volume Integrals, 



§ 23. We have found for the ellipsoid that 



If 
\\\ 






Ydxdydz:=4:7rahc\^^-\--—-\' . . . jV 

where V is any function of x, y, % whose differential coefficients are finite at the 
origin. It is obvious, therefore, that these two results lead to an infinite variety of 
definite integrals. For example, let us use the second result in finding the value of 
IWx^^y'^'^dxdydz, which is an integral evaluated by Lagrange in determining the 
potential due to an ellipsoid (Todhunter's ' History of the Theory of Attractions/ 
vol. ii, p. 158). 

Putting i=Z4"'?^+^? "we find at once 

[[[ 4:7r(ibc V a^W'^^'C^^^ (P^ (P^ d?^ 

- .^....j.....,<>.. M-5 ■ • . (2?-l)1.3.5 . . ■ (2m-l)1.3.5 . . . {2n -l) 
~ 1.3.5 . . . {2l + 2m, + 2n + '3) 

If we suppose the function V such that it satisfies the equation 

then the surface volume integrals are the simplest possible in regard to results, but it 
will depend upon some other condition attached to V what class of function we shall 
have succeeded in integrating. For example, if V also satisfies the condition 

x-z-+y - + % 7 = B V 

dx '^ dy dz 

we should be led to a class of results similar to those obtained in the case of spherical 
harmonics for the sphere, and, in fact, derivable from these by the theory of correspond- 
ing points. As another example, let us suppose that V, besides satisfying equation 
(42), also satisfies Laplace's equation ; then one solution for V will be of the form 



'F{^h^-c^x+^c^-a^y+^a^-})h) 
We thus get 



Fdxdydz= -^^ F(0) 

MDCCCLXXIX. 3 G 



• (43) 
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As a limiting case, let the ellipsoid wear down to the focal ellipse, and for the sake of 
brevity let the semi-axes of that ellipse be denoted by/* and g. We find 



"^{g^+ify) 






r 

f 9' 



dxdy=^ 2iTfg'F{0) 



s. 



V 



x^ y^ 



2TTfg- 



1 ~j,~^ ^{g^+jfy)d^dy= -^f(o) 



3 



J 



• • • • t 



(44) 



These results may be further simplified if the ellipse becomes a circle. 



The Potential due to an Ellipse of uniform density. 



24, The intoOTal 



e 



O.X + ^y, 



'pds 



taken round the perimeter of an ellipse of semi-axes a, h, is derivable from the corres- 
ponding case of the circle of radius R by the theory of corresponding points, and the 



result is 



QO 1 



« • 



• * • • 



(45) 



If we write herein ad, hO, and recollect that j9§^=:^Sp^ then 



Q^x-^^ygxdy 



taken over the area of the ellipse is equal to 

e'-'^-^^y'^dsd 
the limits of 6 being from to 1. The result of integrating this is 






f » * ft »■ «■ 



• • (^6) 



Similar reasoning would give tis 



yii^'t.y''''dxdy 



a^ ' 52 



Confining our attention^ however^ to the result just found (46), we see that the 
integration over the area of the ellipse of any function V of x^ y^ having finite differen- 
tial coefficients at the origin^ leads to 
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^s!!.^. jail I 

where, after the differentiations are performed, x, y, z are to be put equal to zero. 

If we suppose the ellipse coated with matter of unit density we thus s^et for the 
potential at any outside point (/, g, A), sufficiently distant, 

(Q^ — - Jl. J3 — - j 
dx dyj 1 



x/'{f-xy + (g -yf + J>? 22'+i i\i+l\ r 

/yZ j , fl'i — — j 

.. i^A df da 1 



22^+1 ^!^+l! s/p+g^-Vh^ 

» f~~-l)''2^! 1 / 
= 27ra6|-^g.^^^j^._^'^j 7«(^^'^^+ • • • +B3,(2^, 2cr)+ . . . ) . . (48) 

Where A, . . . B^^ . . . are the same as in § 21, provided we write /'^=a^-- W and g^=:W 
in the expansion of that article. We must have r>b and >'\/a^— 6^ if the above 
expansion is convergent. 

If a=h or/=0 the expansion just found reduces to the well known expansion for 
a circular plate of uniform thickness, given in Thomson and Tait, p. 406. 



(Added September 23, 1879.) 

§ 25. The series given above for the potential of an ellipse may be thrown into the 
form of an integral. Writing the result for the elliptic ring in the form 



# .,dn^i 



27rab%i (a^— 6^);7^— 6^ 



o 



df dhj r 



2'Hlil 



we see that this series is the part of the expansion not involving powers of k or its 
reciprocal of 

k'l / d . d \ j /, d d\'t 

27rahe ^ ^ ^^^ ^/ e^^^ ^^^ ^r- 
where o^=#— P. That is, the series is equal to 

t2ir 

ah \ e^^ ««^ "^ir "^ ^ J -dxjj. 

O Gr -4 
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Hence writing herein aO, bO, multiplying by — and integrating, we find for the 
potential at (f^ g, h) the expression 

I ^(i5cos^~-.sin^|)r__5)^5)^^ (a) 

Jo 

The expansion of the exponential leads of course to the series (48) given above, and 
it is easy to find thereby the expressions for the action of an elliptic current on a 
magnet placed at any distance from a coil outside a certain boundary. In like manner, 
subject to a similar restriction, we might find, as in the case of two ellipsoids, the 
action of one elliptic current on another. Within the boundary referred to, the 
expressions we have obtained are no longer convergent, and in the case of the elliptic 
current the potential in the neighbourhood of the centre of the circuit must be found 
by an independent process. 

The case of the chcle, however, admits of a simple and complete determination. 

Potential of Electric Currents in Circular Circuits. 

§ 26. Reverting to the integral (a) let us put c=^0; then the potential at (/ g, h) 
due to a circle of unit density is given by 

Y=aA \ e-^'^^ ^«^ ^7a -6>dl9dii 
Jo Jo 

To integrate this in regard to 6 we observe first that 

'1 pit pU 1 



Ih l(? V? 



Let us next consider the integral 



r2'n- r37r pk cos ^ 

^/^cosi//^^ or cos li -d-ii 

Jo Jq ^COS-v/r ^ 



Integrating by parts we find 



^C COS ^ ^^^ ^ 



27r 







J cos Y J ^0^ r 



TT n n . OTT 







The first term becomes infinite when. \fi=:- and when \jj=:-~~; we may replace it by 
sec^" \jjd\jj. We thus find 
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~-—j^~~y-~^~------^j~j )dxh=-- cos xbe^' ^°^ ^#= sin^ tbe^ ^^^ '^(iti; 

. \f^ cos y 7r cos-^ ■x^ F cos^ '^/ '^ J o * Jo 

Combining these results, we find for the potential at (/, g, h) due to a circular plate 
of unit density tb e expression 



a^ 



sin^t//e'^^^«^'^5i™# ......... (b) 



If we move the origin to a point at distance v on the negative side of the axis of z^ 
we find 

Jo ^ 

in which form the expansion of —--^ readily gives Clerk Maxwell's expansion for 

the potential of unit current in the circle ('Electricity and Magnetism/ vol. ii., p. 301 

(6'))- 

The result expressed by (b) may be written 

Hence the potential due to unit current in the circle is 



# 



'27r // 1 







dh x/f^ -\-g^-\-{h + aj cos '^Y 



r27r ^ 1 

aj\ sin ^~ ''T T^TTT'Z n rr^ r\l# 

'2"^ aj cos -^ , 

%/f^-\-g^-\-{h-\-aj COB '^y ^ 

^^ p^ aj cos ^ J , . 

Jov/ +^ +(^ + ^J cosy)-^ \/ 

This result shows that the potential of a circular current is the same as that of an 
imaginary bar in the axis of z joining the points whose distances from the origin are 
aj Biid —aj, the density at a point \fj being -—2 cot \p. Now the integral (c) is true at 
all points whose distances from the origin are greater than a. We can, however, 
determine the corresponding result for points within the radius a by the ordinary 
theory of inversion. We have seen that for points outside the radius a the potential 
is the same as that of a bar joining two imaginary points. If z be the distance of any 
point of this bar from the origin and z' the corresponding inverted distance, we have 
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. '. rf/=z: —a/ sec ^ tan t|r(it|i 



a 



and the corresponding density is — 2 cot ti— 7 

■*■ ^ *^ ^ ty cos y 

Hence the potential at all points inside the radius a is 



2 



f" ^ se^^ t ^^ ....,,,.. (d) 



This result may also be proved from the consideration that the expression for the 
magnetic force is continuous as we pass from one form of potential to the other. 

§ 27. The results (c) and (d) give the mutual potential energy between two circular 
currents placed in any positions whatever, in which the smaller circle is either wholly 
outside or wholly inside the radius of the other. For the sake of simplicity we shall 
suppose that the axes of the two currents intersect. Let the radius of the larger be 
A and of the smaller a, and let ii, v be the distances of their centres from the point of 
intersection, the axes being inclined at an angle cos""^ /x. Then, if unit current cir- 
culates in each, and if the smaller circle lies wholly inside the sphere of radius A 
containing the larger, the mutual potential energy is given by 

. A f'T'' / cos 6 sec ^-vlr T , J . 

J oi Q's/iu + Aj sec -f) + (v + aj cos ^f — 2(it + A; sec 'f) {v + aj cos (j>)fM 

If this series be expanded in harmonics 

B0P0+B1P1+ . . . +B,P,+ . . , 



we have 



B.=4aA "^ 7 TT , : ^, ^ d6d\b 

'' ' (2^ + Ajseci|r)«+i ^ ^ 



JO^ 



The series in the form of zonal harmonics is given by Clerk Maxwell {' Electricity 
and Magnetism,^ vol. ii., p. 303). 

On comparison with his series, since 



J cos (t>{iM+jp cos <l>Ydcf>= -:^^^^ 
when we take account of the above value of B,, we deduce 



V — " 



o(/jb + vj sec '^)'*"^-'- ^ dfM 

true for integral values of n, from 1 upwards. 
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If the smaller circle is wholly outside the sphere of radius A, the value of M given 
above is inapplicable, and we ought to take a double application of the formula (c). 

§ 28. In the important case of two coaxal circular currents, if & be the distance 
between them^ the mutual potential energy is given by 



M=4aA 



TT rTT 



j COS 6 sec^ '^7,7. / V 

-^ ^ ^ dcjydxjj (e) 



oJ 0^ +i(A sec yfr—a cos cj)) 



provided the sphere of radius A encloses the smaller circle. 
Performing the integration in regard to i// we find 

M=47ra /,, , , ^.M 

"With the view of simplifying this let us put a cos (^+&y=-^cos (x+yj), so that 
a cos ^= A cos cc cos 3^* and — hj=-A siiix sin yj. Thus x and y are not independent, 
and by the substitution proposed it can readily be shown that M takes the form 



SttA 



3 



cos^ X sin X 



\/A^ cos^ x — (a^-hA^-\-b'^) cos^ x-\-a'^ 



dx 



the limits of integration being determined as follows : 

When (^=2^ ^=^^ and when ^=0, A^cos^cc — (a^+A^+&^) cos^ x-\-a^=0. 
The roots of the last equation are 

p^= cos^ 0^ = 2A? 

^3_ ,,,3 fi _ A^ + ^^ + ^^ - y (A^ + 6^^ + ^^) - 46^W 

Since _p is greater than unity we must take q as the value of cos ^ when (f)=-0. 
Hence the value of M is 



»7r 



. 1 2 cos^ X sin X , 

cos~"i q"^^^" ^^^^^y^^" ^^^^^) 



a 6 




a/p^ — q^ cos^ 




^ 



n 



^^M%)-%)\ (f) 



where F and E are complete elliptic integrals whose modulus is -. 
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If h be small and if A — a=:c be also small, then approximately 

-^^ and .2= 1 - V "A^ 
and M becomes 

4.A{log8^^-2 

This result is given by Clerk Maxwell. 

§ 29. If the circle of radius a be outside the sphere of radius A passing through the 
larger circuit, the value of M is 

. . r r cos cf) COS -i/r 

JoJo&+i(Acos^ + (^cos'v|r) ^ ^ '^^^ 

The forms (e) and (g) obtained for M are interesting from their simplicity, and would 
seem to be useful in calculations connected with induction coils. 



